Scaling state of dry two-dimensional froths: 
universal angle-deviations and structure 
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We characterize the late-time scaling state of dry, coarsening, two-dimensional froths using a 
detailed, force-based vertex model. We find that the slow evolution of bubbles leads to systematic 
deviations from f 20° angles at three-fold vertices in the froth, with an amplitude proportional to 
the vertex speed, v ~ i _1//2 , but with a side-number dependence that is independent of time. We 
also find that a significant number of Tl side-switching processes occur for macroscopic bubbles in 
the scaling state, though most bubble annihilations involve four-sided bubbles at microscopic scales. 
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I. INTRODUCTION 

Froth systems have a distinctive, close-packed cellu- 
lar structure of bubbles that are bounded by discrete 
smoothly curving walls (films) that intersect in point- 
like Plateau borders (vertices) . Random soap froths con- 
tinually coarsen due to gas diffusion between individ- 
ual bubbles 0, 0, 0. These dynamic froths exhibit a 
characteristic length-scale that grows as a power-law in 
time, L ~ t 1 / 2 , as well as a time- independent "scal- 
ing" structure. The scaling state, or regime, of froths 
is observed experimentally in both two-dimensional (2d) 
and three-dimensional froths 0. Our understanding 
of the scaling state of 2d soap froths is largely built on 
the foundation of experimental [y, 1?J and computational 
H H E3 [H [H E3 studies close of 2d froths, in the 
dry-froth limit with Plateau borders of vanishing size. 

In the dry- froth limit, and under the approximation 
that bubble-wall curvatures and intersection-angles ex- 
actly satisfy local force- balance, individual bubble areas 
evolve with von Neumann's law, 



dA n /dt = iraD{n - 6)/3, 



(1) 



where A n is the area of a bubble with n sides, a is the sur- 
face tension of the films, and D is the coefficient of diffu- 
sion of gas between adjacent bubbles 0. This slow, con- 
tinuous evolution is supplemented with dynamical rules 
for fast, topological rearrangements: both "T2„" pro- 
cesses of bubble annihilation for n-sided bubbles, and 
"Tl" side-swapping processes between adjacent bubbles 
(see Fig. n below). 

Experimental studies of undrained two-dimensional 
froths exhibit growing Plateau borders and system- 
atic deviations of three-fold vertex angles from the force- 
balance angle of 27r/3 0. Conversely, experiments 
done in the dry-froth limit have not shown measurable 



angle-deviations [6j. The angle deviations reported for 
undrained froths have been previously explained by the 
combination of film curvature and finite Plateau bor- 
ders [Til Hoi Il6j| . However, small angle deviations are 
required, even for dry froths, in order to allow vertices to 
move and accommodate continuous bubble growth. In 
this paper, we use a natural force-based vertex model 
(see Sec. [Hj that allows us to measure angle-deviations 
for dry froths. We find angle-deviations that are quan- 
titatively consistent with experimentally observed angle 
deviations of undrained froths 0. 

Experimentally |17|. four- and five-sided bubbles have 
been observed to annihilate directly through T24 and T2§ 
processes, without any intervening Tl processes, consis- 
tent with theor y |X8 I- However, previous numerical stud- 
ies llllEEElllllHlIIl have not clearly distinguished 
between four- and five-sided bubbles that directly an- 
nihilate at microscopic scales and those that shed sides 
through Tl processes at macroscopic scales before anni- 
hilating. There is a related controversy over whether any 
Tl events occur for macroscopic bubbles in the scaling 
regime 0], or whether all late-time Tl events arc 
associated with microscopic side-shedding during bubble 
annihilation pol l2ll l22j . In the scaling regime of our 
model froth, we confirm the experimental observations 
of direct annihilation of four- and five-sided bubbles at 
microscopic scales. We also observe a significant rate of 
macroscopic Tl events in the scaling state. 

A standard measure of the topological structure in 2d 
froths is the second moment of the side-distribution func- 



tion, fl 2 



where n is the number of walls 
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of an individual bubble. Detailed froth models have re- 
covered H2 values ranging from 1.2-1.6 EJ 13 E3 
with reports of sensitive dependence of \ii on the de- 
tails of T24 and T2s processes [TH The frequently 
quoted experimental value of \ii = 1.4 ± 0.1 comes from 
an undrained froth experiment Measurements of /12 
are bedeviled by long-lasting initial transients. In this 
study, we find /12 = 1-24 ± 0.01 and show that the tran- 
sient after the initial [ii peak can last a decade in time. 
Our result is consistent with the late-time experimental 
results of undrained froths Q , and with numerical work 
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FIG. 1: Illustrations of the fast topological processes involved 
in froth evolution. From top to bottom are a side-switching 
Tl event, and three-, four-, and five-sided bubble deletions 
(T23, T2i, and T2s respectively). From left to right are ex- 
amples of the bubble walls and vertex positions before, during, 
and after the corresponding topological process. The details 
of resolution of T2i and T2$ events will, in general, be deter- 
mined by the forces acting along the bubble walls, as described 
in Sec. El 



on a locally equilibrated dry-froth [l3[ . 

Computational model froths have either been fast but 
simple "topological" models (see, e.g., (3), or slower but 
more detailed "vertex" models with vertex positions and 
film curvatures, H H [H GH 13 E3 With a natural 
force- resolved vertex model (section II), we characterize 
the scaling state of two-dimensional dry froths. Starting 
with up to 10 4 bubbles, we reach the scaling state af- 
ter long transients. In the scaling regime, we determine 
the second-moment of the side-number distribution /l*2 
and the area growth-law amplitude (A)/(tDa) (section 
III. A) . We also characterize the angle deviations at bub- 
ble vertices, and find universal angle deviations associ- 
ated with slow bubble coarsening (section III.B). Finally 
we find a significant number of macroscopic Tl events 
occurring in the scaling regime, and find that most bub- 
bles annihilate as four sided bubbles through T24 events 
(section III.C). 



II. FORCE-BASED VERTEX MODEL 

We implement a fully force-based vertex model with 
natural deterministic dynamics for both the slow evolu- 
tion of bubble areas and for fast topological rearrange- 
ments. The motion of vertices is driven by the vector 
sum of the surface tension of the three films that meet at 



each vertex, leading to vertex velocities of 

v = ia^2 i, (2) 

i 

where the % are unit vectors tangent to the films meeting 
at a vertex and 7 is the vertex mobility. We allow for 
film curvatures, k, to dynamically relax towards their 
steady-state by 

dn/dt = T(AP - no), (3) 

where V is the curvature mobility [2^ and AP is the 
pressure difference between the two bubbles separated 
by the film. Froth coarsening is driven by gas diffusion 
between bubbles with a rate proportional to the diffusion 
constant, D, the full length of the intervening film, and 
the pressure difference across the film: 

dAf/dt = D I APdl, (4) 

where we impose an ideal gas relationship JV = PA, 
where A is the bubble area and P is the bubble pres- 
sure. The continuous evolution of vertices, curvatures, 
and bubble area is implemented with an Euler time- 
discretization with timestep At, supplemented with the 
dynamical rules for topological processes. 

The essence of both Tl and T2 processes is the tempo- 
rary coalescence of two or more three-fold vertices into 
an rt-fold vertex, followed by the dissolution of that n- 
fold vertex. This all happens at microscopic length and 
timescales, determined in experimental systems by the 
Plateau border size. Our approach implements this phe- 
nomenology, though we choose computationally conve- 
nient microscopic length-scales and rely on the univer- 
sality of dry-froth structure to those scales. We enter 
into Tl topological processes when adjacent vertices ap- 
proach within a microscopic distance, r c , of each other, 
and into T2 processes when a bubble area decreases be- 
low lOr;?. A finite r c allows us to use a fixed timestep, At 
|l3j |. For a Tl process, we temporarily replace two partic- 
ipating three-fold vertices with one four-fold vertex at the 
midpoint, while, for T2 n processes, we replace n three- 
fold vertices with one n-fold vertex at the centroid of the 
deleting bubble. To resolve any n-fold vertex, we first 
determine the two adjacent films whose combined force 
upon the vertex is greatest, in the center of force frame. 
A short film of length r c is placed between these two 
films and the other films of the vertex, oriented with the 
combined force, giving an n — 1-fold vertex and a three- 
fold vertex; this process is repeated until only three-fold 
vertices remain. This approach reproduces the natural 
instability of four- and five- fold vertices > an( i follows 
the lowest energy (highest force) channel for resolving 
n- fold vertices. The entire topological process is imple- 
mented between timesteps. Our results are not sensitive 
to r c , as long as it is much smaller than the average bub- 
ble scale, or to the precise details of the placement or 
resolution of the n-fold vertex at scales of r c . 
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Except where otherwise noted, we use a = 1, 7 = 
1, T = 0.5, and D = 0.2. We use a timestep At 

0. 01 and r c = 0.01. Our systems are initialized with 
up to Nq = 10 4 bubbles with a random but periodic 
Voronoi construction. To control for finite-size effects, we 
have investigated different Nq while keeping a fixed initial 
average bubble area and pressure, Aq = 10 3 and Pq = 1, 
respectively. We average at least 50 independent samples 
for Nq = 10 4 , and displayed error bars are statistical |24j . 

It is useful to compare our model to two recent ver- 
tex models, one by Herdtle and Aref 13] and another by 
Chae and Tabor |l2j . Herdtle and Aref build their model 
around an exact implementation of von Neumann's law, 
Eqn. ^ i n the late-time limit and neglect any effects of 
angle-deviations at vertices or of film-curvatures. This 
exact implementation requires computational power of 
0(N 1 - 3 ) for each timestep where N is the number of 
bubbles, as opposed to O(N) for our algorithm. In terms 
of our dynamics, they implement the 7, T — > 00 limit in 
Eqn. (J2J) and © , respectively. They only allow T2 3 bub- 
ble deletions, and resolve shrinking four- and five-sided 
bubbles very finely in time to catch all of the determin- 
istic Tl side-shedding that occurs at microscopic scales 
I Kit . Chae and Tabor also impose curvature equilibrium, 

1. e. r — > 00, but have taken a length-dependent ver- 
tex mobility 7 cx where U are the lengths of the 
films adjoining the vertex. At late times, as typical side- 
lengths grow as I ~ t 1 / 2 , they will approach the 7 = + 
limit — which we expect to be singular. They allow all 
T2 n processes, but impose various side-shedding rules for 
resolving T2 4 and T2 5 processes — and find that /Z2 de- 
pends on which ad hoc scheme is used |12) . Like these two 
vertex models, our model is also effectively in the high 
curvature mobility limit compared to the natural curva- 
ture dynamics 23] and indeed we find that curvatures 
quickly approached the equilibrated limit k — AP/a as 
t — ► 00 (data not shown). We believe that our force- 
based resolution of Tl and T2 n events is equivalent to 
that of Herdtle and Aref ^3 since the same film asym- 
metries that lead to side-shedding at microscopic scales 
in their model will lead to vertex shedding from n-fold 
vertices in our model. 



III. RESULTS 

A. Scaling State 

Two-dimensional soap froths exhibit dynamic scaling 
at late times. In the scaling state, the structural prop- 
erties of the froth are time-independent once appropri- 
ately adjusted for the shrinking number of bubbles N(t) 
and growing average area per bubble A(t). Topologi- 
cal structure such as the distribution function of side- 
number of individual bubbles becomes time-independent 
without scaling. Correlations in cellular systems, such 
as 2d froths, appear to be universal due to the strong 
separation in length- and time-scales between the con- 



tinuous evolution of area and the fast Tl and T2 n pro- 
cesses. Dry 2d froth experiments at various temperatures 
recover universal correlations , as do undrained exper- 
iments that vary the gas phase and hence diffusivity, D 
4J. Our results are consistent with the existence of a 
universal scaling state at late times, though we have not 
systematically explored parameter space. 
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FIG. 2: The average bubble area of a coarsening froth, {A}, 
vs. the natural time variable t = tDa. We see that lin- 
ear growth is obeyed at late times, with a best-fit amplitude 
{A) = 0.94 ±0 - 01 (t-t ) indicated by a solid line This am- 
plitude is expected to universally apply to coarsening froths 
in their scaling state, with initial conditions determining the 
effective time origin to- The same symbol legend applies to 
other figures. 

From von Neumann's law, we can directly see that 
(A) ~ t 0, and this is generally observed [a. I2I l25|. 
Since the froth evolution is driven by von Neumann's 
law (Eqn. ^| , the natural time- variable is 

t = tDa, (5) 

and includes both the diffusivity and the surface tension, 
but not the curvature mobility, r,23j. As shown in Fig.|3 
we recover the expected coarsening law for t > 10 4 , af- 
ter initial transients. We see finite-size effects enter for 
t > 2 x 10 5 for No = 10 4 — roughly when the number of 
remaining bubbles drops below 100. We fit |2J] the am- 
plitude (A)/t = 0.94 ± 0.01 in Fig.0— which we expect 
to be universal. 

We next measure the surface energy of the froth, E, 
normalized by i^ex = Nct\J2\/ZAq, the energy of a reg- 
ular hexagonal froth with the same average bubble size, 
A [H [l|. Our result 0, E/E hcx = 0.938 ± 0.001 
is shown in Fig. [3] and is consistent with the best previ- 
ous measurement (0.945±0.010 0]), though with higher 
precision. 

Simulations of dry froths [1 GH H El HI have given 
values of the second-moment of the distribution of the 
number of sides in a bubble, {12 = (n 2 ) — (n) 2 , that range 
between 1.2 and 1.6. As shown in Fig.QJ we observe [24J 
a constant, asymptotic value of ^2 = 1.24 ± 0.01. While 
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FIG. 3: Film energy of a coarsening froth, E, equal to the 
surface tension a times the total film length, normalized by 

the energy of a uniform hexagonal froth -Ehex = Nay/ 2\/3Ao, 
with the same number of bubbles, N, and average bubble 
area, Ao, plotted vs. t = tDa. The asymptotic value in the 
scaling state is E/E hcx = 0.938 ± 0.001 24] indicated by the 
horizontal line. 
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FIG. 4: The second moment of the side-number distribution, 
fi2 vs t. We extract [24|] an asymptotic value of /12 = 1-24 ± 
0.01, indicated by a horizontal line. Significant transients are 
seen for the larger D data ("+"), including a distinctive peak. 



broadly consistent with previous values, we are closest to 
Herdtle and Aref 0], who found ^2 ~ 1-2. 

While we expect that A/t, E/E\ lc ^, and [ii are univer- 
sal measures of the scaling state, both the initial froth 
configuration and the particular dynamical parameters 
of the model will affect initial transients. By varying D, 
we can simulate the effects of using various gas-phase 
components — much as experimental studies have used 
either helium or air in their froths Q. For larger D, we 
qualitatively reproduce a strong pre-asymptotic [ii peak 
in Fig. 0] that was seen experimentally |4|. The transient 
regime after that /i2 peak lasts approximately one decade 
in time. In the best experimental determination of /X2, 
the entire experiment lasted only one decade after the [ii 



peak |4j, and the latest time data also recovered /Z2 ~ 1.2. 



B. Universal Angle Deviations 

In a force-balanced steady-state, the angles between 
the three adjacent films that form a vertex will all be 
#0 = 2tt/3. Any deviations from 8q will result in a net 
force and motion of the vertex due to Eqn. (J2J. Since 
most bubbles in the froth are continuously growing due 
to Eqn. ifTJl. their vertices must be continuously moving 
- necessitating non-zero angle deviations. We expect 
internal angles typically greater than #0 for growing bub- 
bles (n > 6) and internal angles less than 6q for shrinking 
bubbles (n < 6). 

In the scaling regime, the average bubble area grows 
as (A) ~ t, with a characteristic length L = y/ (A) ~ \ft 
and vertex speed v ~ dL/dt ~ oDt^ 1 / 2 . Small angle- 
deviations will, by Eqn. (J2J, be proportional to vertex 
speeds divided by 7c, so we expect (|A0|) ~ D^ 1 ^ 1 / 2 , 
where A8 = 9 — 9q. As shown in the inset of Fig. [SJ we 
observe this time-dependence in the scaling regime. We 
also collapse the data for different D by scaling by jD^ 1 . 

Different size Plateau borders will lead to different ver- 
tex mobilities, 7, and will in turn require different mag- 
nitude angle deviations for the same vertex speeds. If 
we scale out the magnitude, we should be left with uni- 
versal angle deviations that reflect the universal froth 
structure. Indeed, in our data we find that angle devia- 
tions as a function of the number of bubble sides, scaled 
by the value for n = 5, is time-independent in the scal- 
ing regime for 3 < n < 10 (see Fig. [5j|. Underlining 
this universality, the experimental data from undraincd 
froths |U, indicated by circles in the figure, is consistent 
with our results within error bars. Angle deviations have 
been previously explained as arising from hidden angles 
due the finite size of Plateau borders in combination with 
finite film curvatures (see also 0EH). This expla- 
nation is unsupported by our data. Our angle deviations 
are measured at the vertex center and do not include any 
explicit effects due to Plateau border size. 

Topological events also contribute to locally substan- 
tial angle deviations that are, however, sub-dominant 
in the scaling regime. We can qualitatively understand 
these deviations under the simplifying assumption that 
all angles equal 0q immediately before a topological event. 
As illustrated in Fig.Q] immediately after a Tl event, two 
of the six internal angles of the two separating vertices 
are 7r/3 rather than 2ir/3. Similarly, immediately after a 
T24 event two of the six internal angles about two ver- 
tices are 7r/4, and immediately after a T2§ event two of 
the nine internal angles around three vertices are it/5. 
These large angle deviations will last until a quasi-static 
configuration is reached by moving vertices on order the 
bubble size, L ~ t 1 ^ 2 . As shown in the next subsection, 
the rate of topological events is proportional to the rate of 
bubble annihilation, Rti ~ —dN/dt ~ t~ 2 , so the aver- 
age angle deviation due to topological events is expected 
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FIG. 5: Average angle deviations of n-sided bubbles scaled by 
the n = 5 sided deviations, (A9„) /\{A9s)\, vs. n. The data 
is for t = 27969 (triangles and squares) and t — 158217 ("x" 
and "+"), for D — 0.2 and D = 1, respectively. Also shown 
(circles) is the scaled data from experimental undrained froths 
Q, which illustrates the universality of the scaled angle devi- 
ations. Inset: time dependence of angle deviations of n = 5 
sided bubble vertices, |(A#5}|, scaled by 7D -1 (see text) vs 
t = tDo. The solid line shows the expected t -1 / 2 behavior. 



lengths provides a natural division between microscopic 
and macroscopic, I = y/r c L(i). At late times, I will be 
far from both microscopic and macroscopic scales. We 
count Tl events distinctly if all bubbles involved are big- 
ger than £ 2 , and count T2„ events if the bubble has n 
sides as it shrinks past the scale of £ 2 . Scaling £ by a 
constant factor can significantly change the duration of 
initial transients, but it should not not change results in 
the scaling regime. 
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to be only (|A0|) T ~ t 3 / 2 , which is sub-dominant to the 
effects of continuous bubble evolution. 



C. Topological Processes 

It is known that three-, four-, and five-sided bubbles 
can annihilate directly through T2%, T2 4 , and T2s pro- 
cesses respectively 18], however it has not been clear 
what the relative rates of these processes are. To compli- 
cate the matter, while a three-sided bubble can only an- 
nihilate through a T2% process, a four- or five-sided bub- 
ble can shed sides through Tl processes at macroscopic 
scales and then annihilate through a T2 3 process. There 
have been qualitative reports that there are few if any 
Tl processes that are not associated with side-shedding 
on the way to annihilation [ItI l2l|. and, indeed, simula- 
tions that only allow T2 3 processes are self-consistent and 
yield results comparable to other models 0. However, 
T2 4 and T2§ processes have been reported experimen- 
tally [13 . 

What has been lacking is a clear definition of a Tl 
process as distinct from a T2 process. A natural dis- 
tinction is between macroscopic Tl processes and mi- 
croscopic T2 processes. Tl processes should be counted 
distinctly if they are associated only with macroscopic 
bubbles. Any side-shedding that occurs at microscopic 
scales should simply be counted as the details of the as- 
sociated T2 process. We have a natural microscopic scale 
that corresponds to the Plateau border size, r c , while 
the average bubble size, L ~ i 1 / 2 , provides an appropri- 
ate macroscopic scale. The geometric mean of these two 



FIG. 6: The proportion of T23, T24, and T2 5 processes, for 
three, four, and five-sided bubble annihilation, denoted by tri- 
angles, squares, and pentagons, respectively, vs. t. The filled 
symbols are for D = 0.2 while the open symbols are for D = 1, 
both with No = 10 4 . The number of sides are counted as the 
bubble area approaches microscopic scales, at t = y r c L(t). 
The proportions approach constant values of 9 %, 84 ±2 %, 
and 7 ±2 % respectively, for D = 0.2 [24jl . indicated by hori- 
zontal lines. Significant transients are apparent even at late 
times for D — 1, though the trend is towards agreement with 
D = 0.2. 

We have measured the proportion of bubbles which 
vanish with three, four, and five sides, through the T23, 
T24, and T25 processes illustrated in Fig. The pro- 
portions of these processes have not been reported in the 
literature, despite qualitative reports of T24 and T2s pro- 
cesses in undrained froths 0] • As shown in Fig. these 
proportions are constant in the scaling regime, and have 
the values of 9 ±2 %, 84 ±2 %, and 7 ±2 % for T2 3 , T2 4 , and 
T2 5 processes, respectively. What is surprising is how 
few bubbles annihilate with three-sides. This is qualita- 
tively consistent with how few small three-sided bubbles 
are observed in froths in the scaling state; it also indi- 
cates how important T2 4 and T2s processes are in froth 
dynamics. 

We have also measured the ratio of rates of macro- 
scopic Tl to bubble annihilation T2 processes, Rti/Rt2- 
This ratio has been measured in simulations by Herdtle 
and Aref to be 3/2, but they included microscopic 
side-shedding events so their result provides an upper 
bound to ours. As shown in Fig. [7| we find a ratio of 
Rti/Rt2 ~ 0.53 ± 0.05. This indicates that Tl events 
are significant even in late stages of froth coarsening. 
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FIG. 7: The proportion Rti/Rt2 of macroscopic Tl events 
to microscopic T2 events, where only side- shedding involving 
bubbles larger than £ = y/r c L(t) are counted as Tl events. 
We find a ratio of Rti/Rt2 = 0.53T0.05 in the scaling regime 
|24f . indicated by the horizontal line. Very long transients are 
found for D = 1 ("+"), consistent with the previous figure. 



IV. SUMMARY AND DISCUSSION 

We have accurately characterized the scaling regime of 
a two-dimensional dry froth with natural force-resolved 
topological processes. We find values of [1% — 1.24±0.01, 
(A)/t = 0.94 ± 0.01, and E/E hcx = 0.938 ± 0.001, inde- 
pendent of the gas diffusivity D. Our results for \ii and 
E/Eh cx are consistent with previously measured values 
[T3| though considerably more accurate, while our results 
for (A)/i have not been previously reported. These re- 
sults are consistent with the hypothesis of a universal 
scaling state in two-dimensional froths, independent of 
the microscopic details of the model as long as natural 
topological processes are used [l2j. To affect the struc- 
ture, dynamical features should cither change the asymp- 
totic von Neumann's law or the resolution of topological 
processes. So froth rupture 01 an d shear [23 should 
and do affect structure, but varying the experimental 
temperature, the plateau border size, or the initial con- 
ditions do not. We expect that our basic froth model 
yields universal structure in the scaling state if D, 7, T, 
and a are all non-zero, however a systematic exploration 
has not yet been done. From the precision of our re- 
sults, it now appears feasible to undertake such a study. 
We have shown that deviations from 2ir/3 vertex angles 
have a non-universal magnitude (|A0|) ~ i^ 1 / 2 that is 
a natural consequence of slow bubble growth and finite 
vertex mobility, 7. When scaled by the angle deviations 
for n = 5 bubbles, the scaled angle deviations are time- 
independent and agree quantitatively with experiment 
Q — indicating their universality. This universality re- 
flects the universality of the froth structure, as well as 
the robustness of the topological processes to the small 
angle-deviations that remain in the scaling-state. While 
the angle deviations appear to reflect rather than affect 



the scaling structure of the froth, this may not be the 
case in sheared froths [23. Angle deviations must be 
quite large to affect the resolution of Tl processes, but 
even small deviations can affect the resolution of T2^ and 
T25 processes. We therefore expect that angle devia- 
tions will significantly affect structure in rapidly sheared 
coarsening froths. Furthermore, while the angle devia- 
tions following a topological process are sub-dominant in 
coarsening froths, they may be relevant to sheared froths 
since in those systems topological processes are not well- 
separated in time at late times. 

We have also characterized the topological processes 
occurring in the scaling state. We make a clear distinc- 
tion between T2 events occurring at microscopic scales, 
on the order of r c , and Tl events occurring at macro- 
scopic scales, on the order of yj (A). We find that most 
annihilating bubbles have four-sides (84 ±2 %) while rela- 
tively few have three (9 ±2 %) or five sides (7 ±2 %). We 
also characterize the ratio of side-switching Tl events to 
bubble-annihilating T2 events and find a significant num- 
ber of Tl events for macroscopic bubbles in the scaling 
state, Rti/Rt2 — 0.53±0.05. This is a significant result, 
since many topological froth models neglect macroscopic 
Tl processes. This simplification appears to be unjus- 
tified. We would expect small but significant effects in 
the scaling structure due to these Tl processes. It will 
be interesting to block macroscopic Tl events within our 
dynamics, and compare the resulting froth structure with 
our results. 

Our agreement with the results of Herdtle and Aref 
[l3j| . indicates that our natural force-based resolution 
of T2 processes is consistent with their deterministic 
side shedding within T2 processes at microscopic scales. 
This also indicates that their large vertex mobility limit, 
7 — * 00, recovers the universal scaling state of the coars- 
ening froth. Our agreement with the undrained froth ex- 
periments of Stavans and Glazier Q , including the scaled 
angle deviations, indicates that the small corrections to 
von Neumann's law due to Plateau border broadening 
does not appear to change the structure of the scaling 
state significantly. 

Our force-based vertex model is simple, versatile, fast, 
and accurate. It should be straightforward to study 
froths with Nq up to 10 6 in the near future, including the 
effects of shear. This will allow a more detailed investi- 
gation of universal scaling-state structure in soap-froths. 
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